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First-principles investigation of uranium monochalcogenides 
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We present first-principles investigation of the electronic structure and magnetic properties of uranium 
monochalcogenides: US, USe, UTe. The calculations were performed by using recently developed LDA+U+SO 
method in which both Coulomb and spin-orbit interactions have been taken into account in rotationally in- 
variant form. We discuss the problem of choice of the Coulomb interaction value. The calculated [111] easy 
axes agree with those experimentally observed. The electronic configuration 5/ 3 was found for all uranium 
compounds under investigation. 
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Local Spin Density Approximation (LSDA) Ij based 
on the density functional theory is a widely used method 
for electronic and magnetic structure investigations of 
modern materials. Despite the success in description of 
wide band materials, LSDA fails when applied to de- 
scribe transition metal or rare earth metal compounds. 
For example, it gives a metallic ground state in case 
of 3d insulators (such as CuO, CoO, FeO) or underesti- 
mates an energy gap and local magnetic moments values 
for NiO [2] . Another important drawback of LSDA is an 
underestimation of the orbital moment value L. LSDA 
calculations with the spin-orbit (SO) coupling taken 
into account yield the orbital moment value about two 
times smaller then its experimentally observed counter- 
part [21 13] . While L is small for transition metal com- 
pounds and one can neglect it, this is not the case for 
4/- and 5 /-metal systems where the value of the orbital 
moment is larger than the spin one. 

The problem of the underestimation of L comes 
from the orbital independent nature of the LSDA po- 
tential and can be solved by using LDA+U+SO ap- 
proach, where spin and orbital dependent on-site po- 
tential is provided |10j . It results in orbital polarization 
increasing and, hence, it increases the value of the or- 
bital moment and the magnetic anisotropy energy 
The results of previous theoretical [T5] and experimen- 
tal [7] studies demonstrate that the orbital polarization 
mechanism plays a crucial role in 5/ actinide systems, 
such as US, UTe and USe. 
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In this paper we report the results of first-principle 
LDA+U+SO investigations of uranium monochalco- 
genides: US, USe and UTe. Choice of the screened 
Coulomb interaction is discussed and the obtained re- 
sults are compared with the LDA+SO and experimen- 
tally observed one. 

The investigated uranium monochalcogenides have 
the NaCl-type crystal structure and are ferromagnets 
with Curie temperature 178 K, 160 K and 102 K for 
US, USe, UTe respectively [S]. Despite of simple and 
high symmetric crystal structure (FmSm), the easy 
axes are arranged along the diagonal of the cubic cell 
[111] [51 El [S]. In our calculations we use the tight- 
binding linear muffin-tin orbital approach in atomic 
sphere approximation (Stuttgart LMT047 code) [12] 
with conventional local-density approximation and take 
into account the on-site Coulomb and spin-orbit in- 
teractions (LDA+U+SO) [W. The LMTO basis set 
contains the following states: U(7s,6p,6d), S(3s,3p,3d), 
Se(4s,4p,4d) and Te(5s,5p,5d,5/). The Brillouin zone 
integration has been performed on the 8x8x8 grid. 

The value of screened Coulomb interaction, U, 
and Hund's exchange, Jh, are cornerstones of the 
LDA+U+SO method. While determination of former is 
a complicated issue and depends strongly from screen- 
ing in solids, the later is hardly changed from its atomic 
(ionic) value. The constrain LDA method described in 
Ref. [14] gives the value Jh = 0.48 eV for all compounds 
and it is in good agreement with estimated in Ref. [TU] . 
This value of Hund's exchange will be used through the 
paper. 
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Total energy vs. cell volume for different U values ob- 
tained in LDA+U+SO calculations for US. 



It is well-known that the U value depends strongly 
from number of screening channels taken into account. 
For a free ion the Coulomb parameter, U, is about 
20 eV, while its value in the solid varies from 4 eV to 
10 eV for 3c?-metals. In actinides, the value of the on-site 
Coulomb interaction becomes smaller since 5/ states are 
more expanded in real space than 3d one. Several fea- 
tures that can be calculated using band structure meth- 
ods, e.g. magnitude of local magnetic moment, density 
of states at the Fermi energy for metals or energy gap for 
insulator depends strongly on the Coulomb parameter. 
Therefore, the first important task of investigation is to 
define U in reliable way. In this paper we demonstrate 
different approaches to calculate the value of screened 
Coulomb repulsion. 

One way is to adapt the Coulomb interaction value 
to have certain calculated and experimentally measured 
physical quantities coinciding. For example, T. Shishi- 
dou et al. [15] have shown using a Hartree-Fock approx- 
imation with U as a free parameter that calculated and 
experimental values of the local magnetic moments and 
the direction of the easy axes agree well for [7=0.76 eV. 
In this work we fit an equilibrium volume of US by vary- 
ing [/-parameter. The result is shown on the Fig. [TJ The 
calculated equilibrium volume is close to the experimen- 
tal one at the value of [7=1.25 eV and it will be used 
for the uranium compounds under investigation. 

Another way to calculate the value of screened 
Coulomb interaction from first principles is the super- 
cell procedure [14] that takes into account s-, p- and d— 
screening channels. In the framework of this procedure, 
the calculated Coulomb repulsion is equal to 3.6 eV 
for US. This value is larger than reported by Shishi- 
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LDA+SO partial DOS for uranium monochalcogenides. 
j=5/2 is represented by dashed blue line, j=7/2 is rep- 
resented by solid red line and its sum is a dotted black 
curve. 



dou [15j and have found in this work. To understand 
this discrepancy we will analyze the LDA+SO densities 
of states presented in Fig. [2j One can clearly see that 
,7=5/2 and j=7/2 subbands are well separated due to a 
strong spin-orbit interaction (the value of spin-orbit cou- 
pling A=?? eV) and, hence, for the investigated uranium 
monochalcogenides one more screening channel should 
be taken into account. The additional screening chan- 
nel of conducting j=7/2 states leads to much smaller U 
values which are 0.82 eV, 0.98 eV and 0.75 eV for US, 
USe and UTe, respectively. These values are closer to 
one used in Ref. [15] . 

We have performed two series of calculation i) with 
the value of Coulomb interaction obtained by fit of equi- 
librium volume, [7=1.25 eV, for all three compounds 
and ii) with three different [7 values calculated in super- 
cell procedure which takes into account 7/2-5/2 screen- 
ing. 

We emphasize here that in all our calculations we 
have started with initial directions of spin and orbital 
moments coincided with the crystallographic c-axes. 
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LDA+U+SO partial DOS for US. Four most occu- 
pied orbitals shown only. Density of states in Russell- 
Saunders (LS) representation is shown by dashed black 
line. Red solid line shows PDOS in jj coupling scheme. 

The resulting self-consisted S and L vectors were found 
to be anti-parallel to each other and aligned along [111] 
direction. It coincides with easy axes experimentally 
observed for all compounds [HI El [H] ■ 

In order to analyze different contributions of the 5/ 
states to the resulting densities of states we have cho- 
sen a new local coordinate system with z-axis directed 
along total moment J. The partial density of states for 
4 orbitals with maximal occupation numbers in both 
jj- and LS'-representations are shown in the Fig. [3J 
Based on the obtained densities of states we can con- 
clude that there are three occupied orbitals and one 
partially occupied orbital which forms the peak on the 
Fermi level. The similarity of PDOS shape in both jj- 
and LS'-schemes confirms the applicability of intermedi- 
ate coupling scheme. Another argument to confirm the 
intermediate character of coupling is an occupation ma- 
trix which has a large off-diagonal elements in either jj 
or LS coupling schemes. The magnitude of the largest 
off-diagonal clement in both schemes is just ?? per cents 
smaller than diagonal one. This result can be explained 
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Comparison of LDA+U+SO calculation results ob- 
tained for [7=1.25 eV and experimental data 
(g LS =0.7S, gjj =QM). 
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Comparison of LDA+U+SO calculation results 
(Coulomb interaction values are 0.82 eV, 0.98 eV 
and 0.75 eV for US, USe and UTe, respectively) and 
experimental data (gLS=0.73, gjj — 0.86). 

as a strong competition of a strong spin-orbit interac- 
tion and intra-atomic exchange interaction. Indeed, the 
values of spin-orbit coupling, A=?? eV, is comparable 
with the Hund's exchange, Jy=0.48 eV. Therefore, in 
the situation where the Hund's exchange, voting for the 
Russell- Saunders representation, cannot beat the spin- 
orbit interaction, favoring jj representation, the inter- 
mediate coupling scheme is more desirable for all com- 
pounds under consideration. 

Experimentally observed values of the 5/ magnetic 
moments are much smaller than expected for a free ion 
with / 3 configuration (see Table [2]). We have calcu- 
lated the effective magnetic moments for US, USe, UTe 
using the following expression: A^// = 9^bJ{J + 1) 
where Lande factors for f 3 configuration are <?ls=0.73 
and gjj— 0.86 for LS and jj-couplings, respectively. The 
comparison of the calculated and experimental magnetic 
moments is presented in Tables and El Theoretical 
moments seem to be overestimated but howbeit agree 
reasonably well with the experimental data. This over- 
estimation is a result of the large on-site Coulomb in- 
teraction, [7=1.25 eV. The results of second series of 
LDA+U+SO calculation with [7=0.82 eV, 0.98 eV and 
0.75 eV for US, USe and UTe, respectively, are pre- 
sented in Table [H One can clearly see that there is a 
good agreement between theoretical and experimental 
values. 

To conclude, we have performed first-principles in- 
vestigations of uranium monochalcogenides: US, USe 
and UTe. The choice of Coulomb parameter [7 has been 
discussed. We have demonstrated that supercell proce- 
dure with additional screening channels produces the 
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value of screened Coulomb interaction which gives the 
values of effective magnetic moments in better agree- 
ment with the experimental data. It can be traced to 
the fact that this is an electronic degrees of freedom 
based method while fitting of U to have the calculated 
and experimental volumes coinciding should better de- 
scribe phonon properties of solids. We have analyzed 
LDA+U+SO results in two different coupling basis of 
LS and jj types. It was found that for the studied ura- 
nium compounds the intermediate coupling scheme is 
more preferable. Based on the obtained occupation ma- 
trices and partial densities of states we have concluded 
about 5/ 3 electronic configuration for the uranium com- 
pounds. 
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